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Theoretical description of anisotropic systems, such as layered superconductors and coupled spin 
chains, is often a challenge due to the different natures of interactions along different directions. 
As a model of such a system, we present an analytical study of d-dimensional "nodes" arranged as 
the vertices of a Bethe lattice, where each node has nonzero spatial dimension and is described by 
an O(N) quantum rotor model, and there is hopping between neighboring nodes. In the limit of 
large connectivity on the Bethe lattice, the hopping can be treated by constructing a self-consistent 
effective action for a single node. This procedure is akin to dynamical mean field theory, but 
generalized so that spatial as well as quantum fluctuations are taken into account on each node. 
The quantum phase transition is studied using this effective action for both infinite and finite N. 
The importance of the Perron-Frobenius uniform mode on the Bethe lattice is discussed, and its 
elimination via an "infinite range hopping" term shifts the transition, leading to nontrivial critical 
behavior. We calculate critical exponents and find that the internode hopping reduces the upper 
and lower critical dimensions each by one, indicating that-at least for the purposes of quantum 
criticality-the large number of internode couplings is similar to adding a single extra dimension to 
the theory describing a single node. 



I. INTRODUCTION 

It is often the case that the relevent degrees of free- 
dom in a condensed matter system are best understood 
via an anisotropic description, where interactions may 
be strong in some subset of spatial dimensions, and weak 
or irrelevent in others. The most obvious examples in- 
clude the high temperature superconductors, where bulk 
theories are often constructed by studying an individual 
copper-oxygen plane and completely neglecting interac- 
tions between such planes. In the recently discovered 
iron pnictide and chalcogenide superconductors^, such 
an extreme decopuling makes less sense. In this case 
the interlayer interactions must be taken into account, 
although they can generally be treated on some simpler 
footing than those within a plane. Other examples in- 
clude the study of dimensional reduction in organic con- 
ductors, which can be effectively described as coupled 
one-dimensional chainsP While there are a substantial 
number of theoretical tools at one's disposal for investi- 
gating the properties of these chains or layers, treating 
the hopping between them simultaneously in such a way 
that analytical progress can be made is often a challenge. 
Some insight can be gleaned from dynamical mean field 
theory (DMFT), where the hopping between sites in a 
system is treated in a self-consistent fashion, resulting 
in a theory of decoupled sites which retain quantum — 
but not spatial — fluctuations P This approach is known 
to become exact in the limit of large connectivity, i.e. 
for each site coupled to a large number of nearest neigh- 
bors. In general the self-consistent equations of DMFT 
are complicated and often must be solved numerically. A 
simplification occurs, however, when the sites reside on a 
Bethe lattice, in which case the self-consistent equations 



can be solved in closed form and analytical progress can 
be made more readily. 

Since its introduction in 1935, the Bethe lattice^ has 
played a distinguished role in statistical mechanics and 
many body physics. From a theoretical standpoint, the 
tree-like structure of the Bethe lattice, which inhibits 
loops and ensures that there is a unique path between 
any two nodes, often leads to tremendous simplifi cati ons 
and in some cases makes exact solutions possible!^ Al- 
though there is of course no physical realization of such 
a lattice in Nature, it can sometimes serve as a useful ap- 
proximation to real world systems in certain situations. 
Under some circumstances the behavior of a Bethe lat- 
tice may approximate that of a physically realizable lat- 
tice having the same connectivity. One classic example is 
the low-temperature expansion of the Ising model, where 
one finds that the results at the first two orders depend 
only on the number of sites and the number of nearest 
neighbors for each site, so that at this order the results 
are identical for a (hyper) cubic lattice in d dimensions 
and for a Bethe lattice with connectivity z = 2d)r* Con- 
versely, it is sometimes found that the Bethe lattice ef- 
fectively behaves as a one-dimensional system regardless 
of its connectivity. Indeed, one of the main results of the 
work presented here is that the upper and lower criti- 
cal dimensions for a strongly coupled critical theory on 
a Bethe lattice at large connectivity are each shifted by 
one relative to the theory with no Bethe lattice coupling, 
so that the Bethe lattice structure effectively adds one 
extra dimension. Further evidence for the validity of this 
approach comes from the recent study of layered iron- 
pnictide superconductors, where quantitative agreement 
with experimental results was obtained by approximat- 
ing a superconductor — which can be thought of as a one- 
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dimensional stack of two-dimensional layers — as a Bethe 
lattice of layers in the large-connectivity limit P While the 
underlying reason for the similar critical behavior of one- 
dimensional systems and Bethe lattice models appears to 
hinge on the fact that neither system has loops, a better 
understanding of this correspondence is clearly desirable 
and is a motivation for the present study. In general, one 
must always keep in mind the limitations of the Bethe 
lattice as a model for physical lattices and exercise cau- 
tion in relating results between these two systems. 

In this work we consider Bethe lattice systems in which 
the nodes themselves are of finite spatial extent, as in 
the layered superconductor example mentioned above. 
Throughout this paper, the term "node" shall refer to 
a finite- dimensional object composed of d spatial dimen- 
sions. For example, such a model might describe a set 
of parallel, one-dimensional (d — 1) spin chains that 
are coupled to each other as shown in Figure n\ a), or 
a layered system in which each two-dimensional (d — 2) 
layer is coupled to neighboring layers above and below. 
Such an approach has been considered previously in the 
context of coupled Luttinger chains^!, but the empha- 
sis in those studies was on the application of numerical 
techniques from DMFT to obtain information about the 
phase diagrams of the systems under consideration. As 
we mentioned above, the method has also been applied 
to coupled superconducting layers, where the resulting 
effective theory in two spatial dimensions was used to 
obtain expressions for the fluctuation conductivity and 
magnetization.^ Here we study quantum criticality in 
such systems under more general circumstances and us- 
ing an entirely analytical approach. It is worth stressing 
that the approach presented here is distinct from other 
DMFT extensions such as "cluster" DMFT, in which im- 
purities having some nonzero but finite spatial extent are 
embedded into a lattice. In contrast, the "impurities" in 
the present study are of infinite spatial extent. 

We focus here on the example of the quantum rotor 
model because of its relative simplicity, although the 
main results of this work can be applied to a wide va- 
riety of quantum and classical models. In particular, 
many of the results below can be applied directly to the 
<j> model with large N , which is known to share many 
of the same critical properties as the rotor modelP^ The 
quantum rotor model is known to emerge as a low-energy 
effective theory for certain Heisenberg antiferromagnets 
(N = 3), as well as the Bose-Hubbard model (N = 2). El 
Constraining our discussion to the quantum rotor model 
thus allows us to address basic questions about the effects 
of the internode coupling on critical exponents and the 
upper and lower critical dimensions in a relatively simple 
and general setting. In addition, some mean field results 
for nodes of finite dimension can be directly compared 
to those found for similar models with d = using the 
spectral dimension of the Bethe latticeP^^] 

The remainder of this paper is organized as follows. In 
Section II we introduce our model of coupled quantum 
rotor nodes and self-consistently derive an effective ac- 



tion for a single node in the limit of large connectivity. 
The appropriate limits of this effective action for describ- 
ing the ordered and disordered phases are then discussed. 
In Section III we study the critical point in the large N 
limit and determine the upper and lower critical dimen- 
sions. In Section IV we consider 1/N corrections and 
derive results for the self-energy in various dimensions. 
Finally, in Section V we discuss our results and suggest 
possibilities for future work in this area. 

II. SELF-CONSISTENT ACTION IN THE 
LARGE z LIMIT 

We begin with the action describing a quantum rotor 
model for an iV-component field <f> at zero temperature on 
a set of "nodes" forming a Bethe lattice (or, equivalcntly, 
a large random graph with fixed connectivity z), where 
each node corresponds to a quantum critical system of 
spacetime dimension d + 1 > 1 : 

w = E(^ + ^ hop ) C 1 ) 

i 

The first term describes a iV-component quantum rotor 
on each node i: 

Si =- f d d+1 xtf(x)G^tf(x) - ±- f d d+1 x\ t (x) 

+ % - J d^xX^UxWix), 

where g is a dimensionless parameter controlling the 
strength of quantum fluctuations, the field A enforces 
the constraint (j)f(x)(f)f(x) = 1 at every point in space, 
and there is an implicit summation over a = 1, . . . , N. 
The strength of the fluctuations of the A-component field 
4>{x) are controlled by g in ^2j. As g is varied, the sys- 
tem is knowrJEIII] to exhibit a quantum critical point at 
g = g c . For g > g c the system is in a disordered phase 
with (<fi) = 0, while for g < g c the system realizes an 
ordered phase in which <f> becomes ordered along some 
direction. Throughout this study we use relativistic no- 
tation to describe the zero temperature quantum system, 
with p ~ (w,p), and where we have set a velocity equal 
to one. Although much of what follows will be true for a 
general bare propagator Go, for concreteness we can take 
= + r o ; w ith r o > 0. (Note that r is shifted 
from zero to a finite value upon including the simple tad- 
pole correction to the propagator. It can be checked that 
it makes no difference whether this correction is included 
before or after the large z limit is employed to obtain the 
self-consistent effective action.) 

The second term in ([I]) describes hopping between a 
given node i and its z neighbors, denoted by j: 

s ho P = -| £ / <t» l xtf(x)Mx). (3) 

The goal is now to recast the internode hopping described 
by Q as a modification to the action pi corresponding 
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(a) 



(b) 



(c) 



FIG. 1. Schematic depiction of the self-consistent calculation 
of the effective action (111. (a) We begin with a set of nodes 
on a Bethe lattice, where each node describes a d-dimensional 
system and is coupled to a large number of nearest neigh- 
bors, (b) It is assumed — to be verified self-consistently later 
on — that further neighbors can be integrated out, giving a 
correction to the action of the nodes shown in red. (c) The 
integration of the remaining nearest neighbors is performed 
explicitly, leaving a single node decoupled from its neighbors 
and described by the effective action (11) . 



to a single node. Such a procedure shall be shown to be 
valid in the large z limit. In order to describe the transi- 
tion to the condensed phase (with ((pf) ^ 0) consistently, 
it is necessary to eliminate the zero mode solution on 
the graph, also known as the uniform Perron-Frobenius 
mode, which is a general spectral feature of Bethe lat- 
tices or, equivalently, large random regular graphs. This 
mode contributes an isolated delta function at low energy 
to the density of states on the graph, and in a bosonic 
theory there is a tendency to condense into this mode.^ 
In order to avoid such a condensation and obtain non- 
trivial critical behavior, we introduce a constraint that 
eliminates this uniform mode from the spectrum of the 
graph: 



X>?(*) = o. 



(4) 



This condition is enforced by introducing the field x(x), 
which acts as a Lagrange multiplier. This constraint can 
also be thought of as an "infinite range hopping" that 
has equal amplitude between every pair of nodes, as can 
be seen by adding a term ~ x 2 to the action, completing 
the square and integrating out the auxiliary field. 
Let 

S[<f )l ,\ l ,x} = S l {<j> l ,X l }+izt f d d+1 x X a (x)tf(x) (5) 



be the action describing a single node. We shall now 
implement the procedure described above in a self- 
consistent fashion by letting j = 1,2, ...,z label the 
nearest neighbors of node i = and assuming that we 
can integrate all further neighbors out from the partition 
function. The procedure is shown schematically in Fig- 
ure [I] We assume — intending to self-consistently verify 
this assumption below — that this has the effect of mod- 
ifying the action for the nodes denoted by j in the large 



z limit, i.e. S[4n,x] S[4>i,x]- The partition function 
describing the i = node is then (neglecting the normal- 
ization) 



Z(0) ~ / VxV^V\ a e^ s ^ oM ^ 



n 



T>(f>jT>\je 



(6) 



The next step is to perform the functional integrals 
over 4>j explicitly, thus obtaining the self-consistent ex- 
pression for S[<t>o,x]- Expanding the last exponential in 
^ yields 



E 

n=0 



3=1 



(7) 



When the part of this expression in square brackets is 
multiplied out, a typical term will look like 



d d+1 x 2 ^{x 2 )^ 2 {x 2 



jd+l 



(8) 



Now consider the limit z — > oo, with z > n. There will 
be terms such as ([8| with no repeated indices (ji ^ j 2 ^= 
■ ■ ■ 7^ in)j terms with one repeated index, etc. In general, 
there will be z n ~ m n\/(n — m)\ terms with m repeated 
indices. (Note that terms with a single index appearing 
three or more times are neglected, as these terms will be 
of lower order in z.) The sum in ^ thus becomes 



j-n r, 

^ (n - to)! 

n=0m=0 v ' 



d d+1 x^{x)(^(x)) 



x / d d+1 x / d d +V^(z)<^(*)^V)>^0O 



(9) 



The second summation in the above expression runs to 
n/2 if n is even, or to n/2 — 1 if n is odd. We have 
also performed the functional integrals over 4>j, yielding 
expectation values defined with respect to S[<fij,x], as 
can be seen from Q. The sums in ([9| can be performed 
explicitly, resulting in the expression 



ztoe zt f(4,)t J ^ ^ _ zt 2 e ztl J J ^)4> J J ^ ^ 



Zt J 4>{4>) 
exp(r[^,x]) • 



(10) 



The integration variables and vector indices in this ex- 
pression have been suppressed for brevity, but they can 
be restored by comparing with ([9]), which features the 
same integral expressions. The node indices have also 
been dropped, which is permissible since the layers are 
indistinguishable and the only remaining field that has 
not been integrated out is cf>Q = ip. Note in particular that 
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r is a functional of both (f> and x, with the latter depen- 
dence arising indirectly through the expectation values 
(</>) and {4><j)). 

We thus obtain the following result for the partition 
function and effective action describing node i = 0: 



Z(0) = / V<j>VxD \ e - s W' X ^ 
S[ct>,X,x\=S[ct>,X,x}-mx}- 



(11) 



Ideally, one would like to calculate (<j)(x)) and (<j)(x)<f)(x')) 
explicitly with respect to (111 in order to obtain a fully 
self-consistent effective action in terms of the fields <j> and 
X- However, the unwieldy form of Y makes this pro- 
hibitively difficult, so we shall have to content ourselves 
with various approximations. There are two especially 
interesting cases to consider, which describe the ordered 
((0) 7^ 0) and disordered ((<j)) = 0) phases, respectively. 
These limits are considered in turn in the following sub- 
sections. 



A. Ordered Phase 

Let us first assume that the system is ordered, which 
means that one component of ip has a nonzero expec- 
tation value. Letting <f> = (a, n), we consider the case 
where (a) ^ 0, and if describes the remaining N — 1 
components. From the form of (10 1, one sees that in 



this case it makes sense to rcscalc the hopping amplitude 
as t' = zto, which is to remain finite as z — > oo and 
to — > 0. The importance of this particular scaling to ob- 
taining a nontrivial description of the condensed phase 
has been previously noted in th e context of bosonic dy- 
namical mean field theory!^^ In this case the terms in 
( 10 1 with coefficients ~ zt\ disappear in the large z limit, 



and the effective action in (11) reduces to 



5 orc #, A, x] = A, X ] ~ t' I d d+1 x (r(x)) <f> a {x) 
= J d^x^aG^a+^G^-^X (12) 
+ l -X (a 2 + ttV) + H [z X ° - (a)] a + tt' X a A, 



Note that X cannot be trivially integrated out of (12 1, 



since the expectation value (a) depends on it. Assuming 
a uniform solution, the mean field equations for a and x 
following from ( 12 ) are, respectively, 



0=(r-t'){a)+it'(x°), 
= t'(a) r * M 



S X ° 



(13) 



where r = r$ + i\\n=oo is the mass term for the field 
a, including the saddle point value of the field A, which 
we determine self-consistently in the next section. These 



equations have the two possible solutions 

f(cr)=0, r>0 
\(<j)=iX°, r = 0. 



(14) 



This result indicates that there is an instability to an 
ordered phase at r = 0. Substituting this mean field 



result back into the action (12), one obtains the usual 



action for the quantum rotor model, which is k nown to 
exhibit a second order phase transition at r = 0) n l 15 l 

Some additional remarks regarding this result are in 
order. First, it can easily be checked that if the constraint 
Q eliminating the zero mode is not enforced, then the 
transition to the ordered phase occurs at r = t'. (This 
can be seen from the mean field equation for <f> following 
from (12) with the term containing x{p) se t to zero.) 
This would be undesirable because, as we shall see below, 
when studying the disordered phase it is natural to scale 
the hopping as \fzt a = t = finite for z — > oo; but in 
this case the transition would occur at r = \fzt, which 
is infinite under this scaling, implying that the system 
is always in the ordered phase when the zero mode is 



not eliminated. Second, ( 14 ) was obtained with the first 



scaling, in which t' is equal to a finite constant. When 
approaching the transition instead from the disordered 
side and employing the scaling t = finite, it shall be 
seen that the transition occurs not at r = 0, but at r c = 
It. This is not inconsistent, however, since the latter 
quantity is zero in the scaling limit that was used to 
obtain (14). We therefore see that the overall effect of 



the elimination of the zero mode is to shift the transition 
to the ordered state from r = oo to r ~ 0(t). 



B. Disordered phase 



Now let us consider another limit of (11), in which 
we assume (<f>) = and scale the hopping according to 
t = \fzto, which remains finite as z — > oo and to 0. In 
this case, the effective action reduces to 

•Sdis^, A,x] = sy>,x,x\ 

~%\ dd+lx J dd+lx '$ a ( x ){$ a ( x )$\ x '))$\ x 'h 

The mean field equation for (j> following from this action 



o = *(x Q (^)> 



x' [5{x - x')r - t 2 (4> a {x)tf(x'))] {(j> b {x)) . 



Because in the disordered phase (4> b {x')) — by assump- 
tion, this equation implies that (x a (x)} = 0. This re- 
sult is reasonable, since it effectively says that — at least 
at mean field level — the constraint eliminating the zero 
mode has no effect in the disordered phase. It is also 
easy to see that the same result will hold in a general 
interacting theory, such as a theory with a 4 coupling. 
Quite generally, one sees that (0) = =>• (x) — 0. 
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III. CRITICAL POINT IN N -> oo LIMIT 

With the results of the preceding section in place, we 
can now proceed to characterize the quantum phase tran- 
sition at g = g c , starting from the disordered phase in the 
large N limit. As we saw in the last section, obtaining 
nontrivial behavior requires that we employ the scaling 
in which y/zto = t takes a finite value. After rescaling 
g — > g/N, this effective action (15) becomes 



5eff[0, A] = Ud d +h 



(17) 



The new propagator is related to the original one by the 
condition 



where 



g-^(p) = G^{ P )-t 2 g{p), 



G(x-x') = 



(18) 



(19) 



is the full propagator, defined with respect to the effective 



action (j 1 T > . In writing (17), we have used the fact that 
S ab . [Note that there is no sum on a in (19).] 



Together, (17)— (19) form a closed set of self-consistent 
equations describing a single node on the Bethe lattice, 
where the effects of hopping between nodes are included 
in the self-consistent propagator (Jo- 
in the N — > oo limit, we have G{p) 
can be solved, yielding 



Go (p), and (18) 



Go = 



P 



) 2 - At 2 



(20) 



Here we have used for the bare propagator G^ 1 (p) = 
p 2 + r, which includes the saddle point contribution 
from the field A. Although at large momenta we have 
Go — > Go, it is apparent from (20) that there is no value 
of r for which the propagator diverges at p — 0. Sim- 
ilar non-divergent behavior has been noted in DMFT 
studiesP The lack of a diverging susceptibility in this 
case can be traced to the elimination of the uniform 
Perron-Frobenius mode on the Bethe lattice. Such a di- 
vergence would correspond to the onset of uniform order 
along a given direction, but such order is forbidden by 
the constraint (|4). Rather than condensing into a state 
with uniform order on every node of the Bethe lattice, 
the system condenses into the next-lowest energy mode 
available, which is nonuniform from node to node and 
is specified by some complicated "wave vector" on the 
Bethe lattice. A parallel can be made here to antifer- 
romagnetic systems on an ordinary, say square, lattice, 
where the uniform (i.e. ferromagnetic) susceptibility re- 
mains finite due to the fact that the system does not 
have an instability toward uniform order, but rather a 
tendency to form nonuniform ordered state character- 
ized by some nonzero wave vector. In the self-consistent 
large- z calculation presented here, the Green's function 
( 20 1 plays in effect the role of the uniform susceptibility, 



and the fact that it does not diverge simply means that 



the lowest energy mode on the graph is nonuniform by 
construction. 

Although there is no divergence in the Green's func- 
tion, one can see that there are gapless excitations at the 
critical point by considering the spectral function. Re- 
turning for the moment to nonrelativistic notation and 
analytically continuing to real frequency, the spectral 
function is given by 



1 



27rf 2 



-Im Go (wj k) 

e [u 2 - (k 2 

X y/(-U 2 + 



(21) 



k 2 



One sees from this expression that there is a continuum 
of excitations, with the energy gap A = y/r — r c . This 
continuum of states can be explained by thinking of the 
node described by the effective action (17) as an impurity 



coupled to an effective medium with which the impurity 
can exchange energy. Such a picture is often employed in 
DMFT, in which case the impurity is a single quantum 
site. 

In order to determine the value of r, which depends on 
the coupling <?, we first integrate out the field <f> from the 
action dT7|), leading to 



S eS {\] 



N 
~2 



tr In g. 



iX 



d d+1 xX 



(22) 



In the N — > oo limit, assuming a constant solution 
i\(x) = i — tq and using the expression (20) for the bare 



propagator, one obtains the saddle point equation 



1 

9 



d d+1 p d 
( 27r )d+i Qr 

d d+1 p 



^Go(pY 



1 



(23) 



(27r) d + 1 y/( p 2 + r f - At 2 ' 



In order to examine the critical properties of the sys- 
tem, we consider the p — > limit of the integrand in 
( 23 1 . Then the integrand clearly diverges as r — > r c = 2t 
from above, with the critical value of g determined by 
the condition 



1 

9c 



d d+ h 



(2ir) d + 1 y/pi + Up 2 



(24) 



Noting that the integrand ~ 1/p as p approaches zero, we 
see that the integral is IR convergent and that the saddle 
point equation has a solution for d > 0. This establishes 
d = as the lower critical dimension of the theory. 

Taking the energy gap A = \/5, where 5 = r — r c , to be 
small, we can investigate how it depends on the distance 
from the critical point at g 
obtain 



g c . From (23) and (24), we 



1 

9c 



d d+i 



P 



(2tt 



id+l 



\Jp A + 2r c p 2 



y/p 4 + 2r c p 2 + 26{r c + p 2 ) + S 2 



(25) 



We wish to determine the dependence of the above inte- 
gral on d for spatial dimensions d = 1,2,3. The integral 
can be calculated analytically and expanded for small S. 
Defining 



A~( 5 - 5c r, 



(26) 



where z and v are critical exponents, we obtain the fol- 
lowing results: 



A, d=l 
i] c ~ ^ A 2 In A, d= 2 
A 2 , d=3. 



(27) 



Due to the Lorentz invariance of our theory, the dynam- 
ical critical exponent is given by z = 1, which remains 
true at all orders in 1/N. The first result therefore cor- 
responds to v = 1 in d = 1. The last line corresponds to 
v = \ in d = 3, which is equal to the usual mean field 
value for a system with no internode hopping. Since one 
obtains mean field critical exponents above two spatial 
dimensions, we thus establish d — 2 as the upper criti- 
cal dimension in our theory. Combining this result with 
the discussion below (24), we see that, compared with 



the standard quantum critical theory without internode 
hopping, the upper and lower critical dimensions have 
both been shifted down by one, with the upper critical 
dimension shifting from d = 3 to 2 and the lower criti- 
cal dimension shifting from d = 1 to 0. This seems to 
indicate that, at least for the purposes of criticality, the 
large number of internode couplings effectively act as one 
extra dimension in the theory. 



IV. CORRECTIONS AT FINITE N 

We next evaluate the leading 1/N contribution to the 
propagator (19). In order to obtain the propagator for 



the field A, source terms can be inserted into the original 
action (17), then integrating out cf> from (17) and dif- 



ferentiating with respect to the source terms yields the 
propagator. This procedure leads to the effective action 



S = 



1 ( d d+1 n 

2 J J^^^iPWi-F) 
%\(P)I1(P)H-P) (28) 



N f d d+1 p 



4 ./ (2tt) 
/ / d d+1 p 



d d+i 



+ - 2 Jj^Jj^^P)x( q - P m- q ). 

The inverse propagator for the A field is proportional to 
the polarization operator: 



n(p) 



d d+1 k 

j^GoiP + Woik). (29) 



We now wish to calculate the full propagator Q (p) with 
respect to the action (28) to 0(1/N), which is shown 



FIG. 2. Diagrams contributing to the propagator Q up to 
0(1/N). The solid line represents Qo, and the dashed line is 
the propagator for A. 



diagrammatically in Figure [2] and related to Go(p) by 
(18). We therefore have 



G(p) = Go(p)-Go(pMp)G (p) 



(30) 



where the self-energy E(p) contains the corrections from 
the second and third diagrams in Figure [2j 



E(p) = E(p) 



n(o)7 (2tt) 



^^0(9)%), (31) 



where 



*W 4 / (^n^ + " -<«"] -< 32 > 

Before explicitly calculating the self-energy, let us de- 
termine the form of the full Green's function to 0(1/N), 
as given in ( [19] ), subject to the relation between G and Go 
from the self-consistency condition (18). Inverting (30) 
gives 



G- 1 (p) = Go 1 ( P ) + np) + -- 



and combining this with ( 18 ) yields 



g- 1 ( P )^G^ 1 (p) + E(p)-t 2 g(p). 



(33) 



(34) 



This can be solved to obtain a self-consistent expression 
for the full propagator: 



G(p) 



p 2 +r + E(p) + y/[p 2 + r + E(p)] 2 -At 2 



,(35) 



which is valid to 0(1/N). As in the N — oo case, we see 
that the Green's function does not diverge for any value 
of r, which again is a consequence of our elimination of 
the uniform mode on the Bethe lattice. 

We now turn our attention to the self-energy, given in 
(31 ) and (32). Since the last term in (31 1 is momentum- 



independent, it simply gives a correction to the value of 
r. Such a shift does not affect the critical properties of 
the system, however, and so we focus on the momentum- 
dependent part of the self-energy, E(p). Furthermore, 
this shift in r can be omitted in the following calculation 
of E(p), since this 0(1/ N) shift would lead to a contri- 
bution at C(l/iV 2 ) in (3|). 

In order to calculate the self-energy, the polarization 
bubble 11(g) must first be evaluated. Using (20), (29) 



7 



becomes 



U(q,r) 



(2t) 



(d+l)/2 



t 2 



d d+1 k 
(27r) d + 1 



k 2 



(k 2 +f) 2 - 1 



(36) 



(k + q) 2 + f~ \] [{k + q) 2 + f\ 2 - 1 



where we have redefined the momenta as k — k\/2t and 
q = q~\/2t, and rescaled r = 2tf. Careful examination 



shows that (32) will be dominated by large-g behavior, 
so we consider II (q) in this limit. Although the integral 
in ( 36 ) cannot be evaluated analytically, we can obtain 



approximate expressions in various dimensions by using 
the expansions 



(k + qf + f- 



\/[(k + q) 2 



- 1 = 





Y^ 

W+i 


-0 


k 2 — r 


( 2fc 2 ^ 


\d+i 


1 — 1 


q 2 — f 



I 

2«3 4 
1 

2 A- 1 



f » fc 2 (37) 
fc 2 > q 2 . 



Here we have neglected terms that integrate to zero and 
replaced (q ■ k) 2 / (q 2 k 2 ) —> l/(rf+l), as appropriate when 
averaging over angular integrals. Careful analysis of the 



integral in ( 
following be 



36j) for both large and small k reveals the 
raviors of H(q) at large q: 



U(q oo) 




(38) 



In obtaining these expressions, we have again used q = 
q/\f2t. Numerical factors have been dropped since our 
only goal thus far has been to obtain the asymptotic be- 
havior of q in the large q limit. By comparing with (32 1, 



it can be seen that £(p) will be UV divergent in spatial 
dimensions d = 1,2, but finite in d = 0. 

With this information in hand, we can obtain a more 



precise result by performing the integral in ( 36 1 numeri- 
cally for several large values of q, then fitting the result 
to the appropriate form from (38). Due to the weak de- 
pendence of the expansion (37) on r, it is found that the 



precise value of r doesn't affect the large q behavior as 
long as r < A 2 . We therefore set f = 1 in the integrand 
and proceed to calculate ( 36 1 numerically. Specializing 
to d = 2 spatial dimensions, the integral is performed for 
many values of q ^> 1, and the results are fitted to the 
expression 

C 

7 

q 



n(?) 



(39) 



where C is a fitting parameter. This is done for several 
values of the upper cutoff A, and the resulting values of C 
are shown in Figure [3] Extrapolating to the A = oo limit, 
the result is seen to be consistent within our numerical 
precision with C = 1/8, the standard result for a d = 2 
system with no internode hopping. 



c 



0.1250 ; 



0.1245 - 



0.1240 - 



0.1235 - 



0.0 



2t/A (10" 



FIG. 3. Values of the coefficient of I1(ij — > oo) from p9[ ) for 
several values of the ultraviolet cutoff A. . Error bars are 
roughly the size of the points in the figure. 



With these results, we can proceed to calculate the 
self-energy. Using ([39]) , we have (still in d = 2) 



S(p) 



At 

N 



d 3 q 2q 
(2tt) 3 ~C 



[(p + q) 2 +f\ 

At f d 3 q f 
CN J (2tt) 3 3^3 
p 2 \n(A/p) 

3ir 2 CN ' 



(q 2 



(40) 



where the expansion (37) was again used in obtaining the 
second line. If C = i from (39), this agrees exactly with 



the standard result for d — 2 with no internode hopping. 
This correspondence is again related to the fact that f 
does not appear at leading order in the expansion of the 
integrand in ( 40 ) . In the standard theory, one would in- 



terpret the coefficient of the momentum-dependent part 
of the self-energy as the anomalous dimension 77, since one 
could express the propagator as G(p) ~ l/p 2 ~ n . One sees 
from the form of (35), however, that, while this interpre- 
tation applies at momenta p 2 3> r, for small momenta 
one cannot define an anomalous dimension due to the 
fact that the Green's function remains finite rather than 
diverging as a power law as p — > 0. 

Turning briefly to the d = 1 case, we can use ( 38 ) and 



follow a similar procedure to obtain the self-energy 



S(p) 



2f r 



In In 



A 
2t 



— In In 



V 
2t 



(41) 



This double log divergence is also found in the standard 
d = 1 case with no internode hopping.^ While it does 
contribute a (very mildly) diverging term to the self en- 
ergy, it does not lead to an anomalous critical exponent 
in the usual sense. 

To summarize, in this section we have derived the self- 



consistent expression ( 35 ) for the full propagator includ- 
ing self-energy corrections. The self-energies in d = 1,2 



dimensions were then computed at leading order in 1 /N 
and were found to have the same form as the standard 
results from the case where there is no internode hop- 
ping. This is especially interesting in light of the results 
of Section III, where it was shown that various quantities 
including critical exponents, the upper and lower critical 
dimensions and the position of the quantum critical point 
all depend on the internode hopping. 



V. DISCUSSION 

We have presented a self-consistent, analytical treat- 
ment of coupled nodes of spatial dimension d > on a 
Bethe lattice. We have derived and studied an effective 



action (11) which recasts the system as an ensemble of 



decoupled nodes, allowing one to investigate the critical 
properties on both sides of the quantum phase transition. 
We found that there is a nontrivial quantum critical point 
once the zero mode has been eliminated and the intern- 
ode hopping is scaled appropriately. Although there are 
gapless excitations at the critical point, there is no di- 
vergence in the susceptibility, which is a consequence of 
the suppression of the zero energy mode. A calculation 
of the critical exponents at the transition indicates that 
the upper and lower critical dimensions are shifted down 
by one, to two and zero spatial dimensions, respectively. 
This manifestation of the large connectivity limit is in 
contrast to the expectation from the spectral dimension 
of the Bethe lattice that would predict d — to be the 
upper critical dimension 12 . It is interesting that, despite 
this change, the self-energy appears not to be affected by 
the internode hoppings. This seems to be due to the fact 
that, while the calculation of v involves only behavior 
at small momenta, the calculation of the self-energy in- 



volves integrals that are ultraviolet divergent, and these 
divergences are unaffected by the internode hoppings as 
long as A ^S> t. 

One can imagine many possible extensions to the work 
presented here. For example, one could investigate how 
these results might be different for internode hopping 
along other, more realistic lattices rather than the Bethe 
lattice. This may require some numerical work to evalu- 
ate the self-consistent equations [which cannot in general 
be expressed in closed form as in (fl8|], but it would be 
interesting to see the extent to which results that we have 
derived for the highly idealized Bethe lattice would carry 
over to physically realizable systems. It would also be de- 
sirable to generalize to finite temperature, and to extend 
our results to fermionic systems, in which case contact 
could be made with the many results from DMFT relat- 
ing to the fermion Hubbard model and related models. 
Investigating the full effective action ( 11 ) numerically for 
finite z may allow one to obtain a more complete de- 
scription of the quantum phase transition. One could 
also numerically study the problem of coupled extended 
nodes on a Bethe lattice (or, equivalently for most pur- 
poses, a large regular random graph) directly for finite 
z and compare the results with those derived from the 
self-consistent, infinite- z theory presented in this work.^ 
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